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Hyperfine spfittings (HFS) are calculated within the Field Correlator Method, taking into account 
relativistic corrections. The HFS in bottomonium and the Bq (q=n,s) mesons are shown to be in full 
agreement with experiment if a universal coupling uhf ~ 0.310 is taken in perturbative spin-spin 
potential. It gives Af(B*) - M{B) = 45.7(3) MeV, M{B*) - M{Bs) = 46.7(3) MeV {nj = 4), 
while in bottomonium AHF{bb) = M(T(9460)) - M{-qt{lS)) = 63.4 MeV for n/ = 4 and 71.1 MeV 
for nf — 5 are obtained; just latter agrees with recent BaBar data. For unobserved excited states 
we predict M{T{2S)) - M{rib{2S)) = 36(2) MeV, M{T{ZS)) - M{ri{3S)) = 28(2) MeV, and also 
M{B*) = 6334(4) MeV, M{Bc{2S)) = 6868(4) MeV, M{B*{2S)) = 6905(4) MeV. The mass 
splittings between D(2 ^5i) - D(2 ^^o), Ds{2^ Si) - Ds{2^ So) are predicted to be ~ 70 MeV, which 
are significantly smaller than in several other studies. 



I. INTRODUCTION 

Spin-spin interaction in mesons has been studied in a large number of theoretical papers [l]-[6], however, up to 
now some characteristic features of this interaction are not fully understood. This statement can be illustrated by 
theoretical failure to explain two experimental facts: rather small -0(3686) — r]c{2S) mass difference: M(V'(3686)) — 
M{ric{2S)) = 49 ± 4 MeV [7]-[9] and, on the contrary, unexpected large HFS in bottomonium, which follows from 
the mass M(r/t(15)) = 9391.1 ± 3.1 MeV of the r]t{lS) meson, recently discovered by the BaBar CoUab. [10]. The 
rjb meson was observed in the radiative decays, T(3S') — > 777f,(lS') and T{2S) — ^rjbi^S) [10], and later confirmed 
by the CLEO CoUaboration, also in the radiative T{3S) — ?> jr]b{lS) decay [11]. The measured value, A//i?(66) — 
M(T(1>S')) — M{rjt{lS)) = 69.9 ±3.1 MeV, is significantly larger than in most theoretical predictions, thus illustrating 
that modern understanding of HF interaction in QCD remains incomplete. 

In perturbative approach a spin-spin potential between heavy quarks contains the factors, like the strong coupling 
and quark masses, which differ in different models and as a result, theoretical predictions for AnFibb) = Af(T(9460)) — 
M{rib{lS)) vary in wide range: 35 — 90 MeV [1]- [6], [12], being in most cases smaller than experimental number. 

On fundamental level spin-spin potential T4s has been recently studied in quenched QCD on large lattice [13], where 
this potential was shown to be compatible with zero at distances r > 0.30 fm (for unknown reason at smaller r it has 
negative sign with a large magnitude). Although the lattice HF potential remains undefined at small r, its behavior 
at larger r is in agreement with widely used Fermi-Breit potential containing (5^(r) [14]. What is important that in 
lattice QCD, as well as in Field Correlator Method (FCM) [15]-[17], a spin-spin potential is described by universal 
functions, expressed via the field correlators. Moreover, in [17] it was shown that nonpcrturbative HF potential can 
give not small contributions to HFS. 

On the other hand in Ref. [1] a smearing procedure for the tJ'^ (r)-function was shown to be very important, giving 
a large Gaussian smearing parameter for heavy mesons, containing a b quark, so that for the Bq {q = n, s, c) mesons 
and bottomonium a smearing occurs at very small distances and for them the use of the Fermi-Breit potential may 
be a good approximation. For the mesons, not containing a 6— quark, a smearing parameter is essentially smaller, 
both for light mesons and for the D,Ds mesons [1]. However, such dependence of a smearing parameter on a quark 
content does not agree with the lattice and FCM representation about a spin-spin potential as a universal one in 
static approximation, where is defined by the field correlators with universal parameters [13], [17]. 

Recently HFS in the Bq mesons and quarkonia have also been calculated in lattice QCD [18]- [22] and their results 
we shall shortly discuss in our paper. Here we study HFS with the use of FCM, where both perturbative and 
nonpcrturbative spin-spin potentials are presented in analytical form [16], [17] and it allows to analyse the role of 
different physical parameters, defining HF structure. However, a comparison of our and lattice results is rather difficult, 
because in lattice calculations perturbative and nonpcrturbative spin-spin effects are not separated and a characteristic 
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value of the strong coupling anp is not discussed. On the contrary, our analysis shows that in bottomonium and 
the Bq (q = n, s, c) mesons HFS can be described within only perturbative approach, since nonperturbative spin-spin 
potential gives a small contribution. 

It is important that in the Bq{q = n,s) mesons and bottomonium a good agreement with experiment is reached 
taking a universal coupling, aHF = 0.310 [23]. We would like to stress that this number is significantly larger that 
that prescribed in pQCD, where as{mb) ^ 0.18 and as{mc) ^ 0.23 are used; just because of such small coupling small 
HFS was predicted in bottomonium in [5], [12]. 

However, for charmonium and the D, Dg mesons their HFS turn out to be by ~ 15% and ~ 30% smaller than in 
experiment, if the same anp = 0.310 is taken. It can occur for two reasons: if for those states nonperturbative HF 
potential gives essential contributions, or higher order corrections are not small, as it takes place in fine structure 
splittings of the IP charmonium multiplet [24]. The situation is diff'erent for the charmonium excited states, for which 
just with a universal coupling, anp ~ 0.31, a good agreement with experimental HFS: Af(?/j(3686)) — M{ric{3637)) = 
49 ± 4 MeV [7] is obtained. Notice that the scale, corresponding to aRpilj) = 0.31 ~ a-silA is rather large, ji ~ 1.7 
GcV. 

In theoretical models two typical choices of aHP are used: 

1. First one, when "a universal" anp is used. For example, in [2] anp = 0.36 was taken from the fit to the mass 
difference, M{J/il!)- M{r]c{lS)) = 117 MeV; then for this choice predicted HFS in bottomonium, M(T(9460))- 
M{r]b{lS)) = 87 MeV, has appeared to be by ~ 25% larger than experimental number. In [3] a smaller universal 
uhp = 0.339 was used in the heavy-light mesons; however, it is difficult to compare our and their results, 
because in [3] a large string tension, a = 0.257 GeV^, was used, while here (as well as in [1]) the conventional 
cr = 0.18 GeV^ is taken. 

2. Second choice is mostly used in pQCD [5], [12], where a scale jj, = rnq depends on a heavy quark mass and 

therefore the value of anpip) — as{rnQ) is essentially smaller. In [12], as well as in the EFG paper [4], just 
due to the choice of aHpijnb) = 0.18 small HFS were obtained in bottomonium (although the w.f. at the origin 
from [12] have provided a precision description of dielectron widths for T{nS) (n = 1, 2, 3) [25]). 

In FCM a spin-spin potential takes into account relativistic corrections and the current masses are used for a light 
quark, m„ ^ 5 MeV (n = u,d), and nis ~ 200 MeV for a s quark (about a choice of nis see [26]), so that the B,D, 
and Bg, Dg mesons can be considered on the same footing as heavy quarkonia and the Be mesons. 

We shall show here that HFS are sensitive to the value of the vector QCD constant Ay(n/), defining a vector part 
of a static potential in coordinate space. In its turn this constant is expressed via Aj^girif) [27], which at present are 
known with a good accuracy only for Uf = 5 and with 10% accuracy for n/ = 3, 4 [7]. To fix Ay(n/) we assume here, 
as well as in [1], that in the one-gluon-exchange potential (OGE) the freezing value of the vector coupling ay(r)(n/) is 
the same for n/ = 3, 4, 5 (it is denoted as acrit)- Due to such an assumption the HFS dependence on n/ is weakening, 
with an exception of the bottomonium ground states. 

We also calculate HFS and the masses of undiscovered yet mesons: r]b{2S), r]b{3S), B*{1S), and the masses of 
Bg{2S),D{2S),Ds{2S). 

The paper is organized as follows. In Section II the spin-spin potential is given in the form, where relativistic 
corrections are taken into account, as it is prescribed in FCM. Also relativistic string Hamiltonian is presented. In 
Section HI the details of the static potential are discussed. In Section IV calculated w.f. at the origin and HFS for 
the Bg mesons and bottomonium are given and their dependence on the number of fiavors is discussed. In Section V 
a choice of the strong coupling for the D,Ds mesons, and charmonium is discussed. Conclusions of our analysis are 
given in Section VI. In Appendix A the conventional formula for the pole mass of a heavy quark is shortly discussed 
and in Appendix B we describe the self-energy contribution to the meson mass, which is important for heavy-light 
mesons. 



II. THE HF POTENTIAL IN THE FIELD CORRELATOR METHOD 



The conventional form of the Fermi-Breit potential [13], 



Vgg{r) = siS2— r , 1 

9 mxrai 

is widely used in heavy quarkonia, as well as in many nonrelativistic models. It contains the constituent quark masses 
mi and m2, which are model-dependent and can differ by 30%, or even larger, in different models, e.g. the mass of 
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a c— quark, rric = 1.48 GeV, was taken in [2], while in the Corneh potential a larger value, rric = 1.84 GeV, was used 
[28]. A constituent mass is supposed to be the same for all nS and nL states. 

In Eq. (1) the strong coupling OHpifJ-) can differ from the QCD strong coupling as{fJ.) (in the MS renormalization 
scheme) due to higher order perturbative corrections. These higher order corrections in one-loop approximation were 
calculated for heavy quarkonia [29]: 



as{mu) 
p(n/) 



(2) 



but remain unknown for heavy-light mesons, containing a light (or a strange) quark. Therefore in general case the 
coupling auF in Eq-(l) should be considered as an effective one. Notice that its value is smaller than a freezing 
constant of the vector coupling ay(r), which defines the OGE potential at large distanses (or at small momenta) (see 
Eq.(17)). In heavy quarkonia with m\ = m2 = tuq the factor p is known [29]: 

and appears to be small: ^ 6 — 8% for n/ = 3, ^ 3 — 4% for n/ = 4, and < 0.1% for Uf = 5; still in some cases these 
corrections can improve an accuracy of calculations. However, since they are not defined for heavy-light mesons, here 
in all cases wc; will consider anp as an effective coupling, which is factually a fitting parameter. 

The important role of relativistic corrections, even for the Be meson, has been underlined in [1], [3], and also in the 
lattice calculations of B* in full QCD [21]. In FCM relativistic corrections are taken into account in two ways: firstly, 
through the kinetic energies of a quark and antiquark, which enter a spin-spin potential [16], [17]: 

Vss{r) = S1S2— — — (5(r). (4) 

y CO1UJ2 



For this potential a HFS is 



A,,inS) = l^^\R^m^ (5) 
y UI1UI2 



where relativistic corrections are taken into account via the averaged kinetic energies coi{nS),u)2{nS): 

uji{nS) = (y^p2 +ml)ns, u}2{nS) = {\J v"^ + ml) nS , (6) 

which are well defined. By definition they depend on the quantum numbers of a given state nS, growing for larger nS 
states. The important point is that in (6) the masses mi,m2 are not arbitrary (or fitting parameters): they are equal 
the pole masses nic. nii, in heav\^ quarkonia, which are now known with an accuracy ^ 70 MeV for a h quark and 
~ 100 MeV for a c quark (see [7] and references therein). In leading order the pole masses ttiq do not depend on a 
number of flavors, while in the order (as(mQ)^) they slightly depend on n/ (as in Eqs.(A.l) and (A. 2) in Appendix A). 
For heavy quarks we take the following pole masses: rric = 1-41 GeV, rrib = 4.79 GeV for n/ = 4 and mi, = 4.82 GeV 
for n/ = 5. 

For a light quark (n = u, d) we use the current mass m„ = 5 MeV and mg = 200 MeV for a s-quark. The mass 

of a s quark is relatively large (close value of mg is used in the Dirac equation in [3]), because the spectra of the 
Dg, Bs mesons are defined at the scale, ^ <1 GeV [26], which is smaller than the conventional scale 2 GeV, for which 
ms{2 GeV) ^ 90 MeV [7]. 

For excited states the kinetic energies (6) increase and therefore HFS, calculated with the HF potential (4), are 
always smaller than those for the Fermi-Breit potential (1) with fixed (constituent) masses. 

Other type of relativistic corrections enter via the w.f. at the origin, which together with the ojj are calculated from 
the relativistic string Hamiltonian (RSH) H^, also derived in FCM [30], 

The variables enter Hq as the kinetic energy operators. However, while a HF interaction (as well as any spin- 
dependent potential) is considered as a perturbation, then in (4), (5) Wj should be changed by the matrix elements 
(m.e.) (6) [16]. 
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^Prom RSH a simple expression follows for a spin-averaged mass M{nS) [31]: 

M{ns) = — + —^ + — + -^+ EnsiuJred) 



(8) 



In (8) an excitation energy Ensii^red) depends on the reduced mass: LOred = ^^^^^ ' bottomonium the mass 

formula (8) does not contain any overall constant, while for heavy-light mesons a negative (with not small magnitude) 
self-energy term, proportional to {q = n, s), should be added [32] (see the expression (A. 3) in Appendix B). 

We use here the Einbein Approximation (EA), when the variables LtJi{nS), the excitation energy Ens{i-^red), and 
the w.f. are calculated from the Hamiltonian (6) and two extremum conditions, =0 (i = 1, 2), which are put 

on the mass M{nS) [31]: 



coi ui2 m, 

— + — H H + 

2 2 2a;i 2u)2 '^^red 



+ VB{r) 



Vns{r) = E{nS)ip„L, 



(9) 



a;,(nS) = m,-2a,,^^^^^^ ir = l,2). (10) 

Before to define the w.f. at the origin and u)i{nS) in next Section we shortly discuss the static potential Vsir) in the 
Eq.(9). 



III. THE STATIC POTENTIAL Vsir) 

In a Hamiltonian approach a choice of a static potential Vsir) is of a special importance; we take it as a sum of 
linear confining term and the OGE -type term: this additivity of a static potential is well established now in analytical 
studies [33] and on lattice [34], [35]: 

Vn{r)=ar+'-^. (11) 
3 r 

For the string tension we use the conventional value, a = 0.18 GeV^, for all mesons (if their sizes are less than ~ 1 
fm. This our choice is in contrast to that in [3], where in the Dirac equation large a = 0.257 GeV^ was used for 
heavy-light mesons. 

The OGE term contains the vector coupling ay (r), taken here in a particular case from the background perturbation 
theory (BPT) and denoted as aB(r) [36], [37]. Two important conditions have to be put on a vector coupling: 

1. i) As in pQCD, it has to possess the asymptotic freedom (AF) property; just due to this property a static 
interaction depends on a number of flavors. Also the AF behavior strongly affects the w.f. at the origin. 

2. ii) The vector coupling freezes at large distances. The property of freezing was widely used in phenomenology 
[l]-[3] and confirmed in lattice calculations of a static potential [34], where a freezing property was assumed at 
rather small quark-antiquark separations, r > 0.2 fm. 

On phenomenological level the freezing phenomenon has been suggested in [38], where in momentum space the 

logarithm In in as(gr^) was changed by In ^ ^t™" ; thus introducing a regulator 4mg. The mass rUg was interpreted 
as an effective gluon mass, although a meaning of mg is not well defined, since in QCD a gluon has no a mass. Later 
a freezing phenomenon was studied in BPT [37], where this regulator was shown to be equal a mass of the lowest 
hybrid excitation (called the background mass), with Mb = 1.0 ± 0.05 GeV [39] for n/ = 4, 5 and a larger value, 
Mb ~ 1.5 GeV for n/ = [36]. As in [1], we shall call a freezing constant a critical one and denote it as acrit- 

Unfortunately, the critical constants, calculated from the static potentials on lattice, are significantly smaller than 
those in phenomenology and BPT. Here in BPT we use rather large aB{crit) = 0.58 — 0.60 for n/ = 4, 5, which are 
close to acrit = 0.60 (for any n/) taken in [1]. On lattice the freezing effect occurs at small distances, r > 0.2 fm, and 
small acrit{lat) ~ 0.30 in full QCD (nj = 3) [35] and acrit{lat) ~ 0.22 in quenched calculations («,/ = 0) [34] were 
obtained. The reasons of these discrepancies arc not established yet. 

While the critical value acrit is fixed, then with the use of Eq.(17) the constant h.B(nf) (for a given n/) can be 
defined. It is important that this constant cannot be considered as a fitting parameter, because it is expressed via 
the QCD constant ^j^-ginf) in the MS renormalization scheme [27] (see the relation (16)). Therefore one can state 
that Ajgg indirectly defines acrit- 
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In the OGE term (11) a veetor eoupling in coordinate space Q!s(r) is defined through the vector couphng aB{(p) 
in the momentum space [25], [36]: 

oo 

2 /■ , sin(gr) , ^ 
"B(r) = - / dq^^ asiq). (12) 



Here the vector couphng asiq^) is taken in two-loop approximation, 

47r / /3iln) 



0'B{q) = ^-—] , (13) 



where the logarithm. 



*. = ^. (14) 



contains the constant AB(n/) defined via the QCD constant h.-j-j^{nf). The relation between them has been established 
in [27]: 

Ai?(n/) = Ajggcxp I I , (15) 



2/?o 

with /3o = 11 — |n/ and ai = ^ — -^n/. From the relation (15) one can see that for a given n/ the constant AB(n/) 
is significantly larger than A-p^: 

A^j^) = 1.3656A(^ (r^/ = 5); 
A^j^) = l.4238Ag^ (n/=4); 

A^^^U 1.4753A(1^ (n/ = 3). (16) 



At present the QCD constant A^^ (for n/ = 5) is known from experimental value of as{Mz) = 0.1182 ± 0.0012 
[7]. Then in two-loop approximation it gives Ap^(two-loop) = 232(12) MeV. The QCD constants h-jjg for n = 3,4 
are extracted from experiments with lower accuracy, ~ 10% [7]. To define them we fix here the freezing value OLcriti 
assuming that they are the same for n/ = 3,4,5. Then from acHt (17) one can calculate all constants h.B{nf) and 
then from (16) to define Aj^{nf). 

Notice that the mass Mb depends on a, being proportional to y/a, and for a = 0.18 GeV^ the value, Mb = 
1.0 ± 0.05 GeV, was extracted from a detailed comparison of the static force in FCM and lattice QCD [39] and also 
the analysis of the bottomonium spectra [40]. Here we use Mb = 0.95 GeV. We also use here two values of A.|^, 
equal 236 MeV and 245 MeV, which give the critical couplings 0.58 and 0.605. 

From (12) it can be easily shown that the critical couplings in momentum and coordinate space coincide, aB{crit) = 
ctB{Q^ =0) = aB{r ^ oo), and it is given by the expression [36]: 

r ■ \ I \ / ^\ 47r / /3i lnto\ 

aB{crit) = aB{r oo) = aB{q = 0) = — — 1 - — , (17) 



with to = tB(9^=0) = /n(^). 



From (16) and A^ = 245(236) McV one obtains that in two-loop approximation a^iMz) = 0.1194(0.1188), which 
agrees within an error with the world average as{Mz) = 0.1184 ± 0.0012 [7]. In Table 1 we summarize the values of 
Ab,Ajys for n/ = 3,4,5. 

The solutions of the coupled equations (9), (10), like the spectra and w.f at the origin, have been checked in 
numerous studies [23]- [26], [39]- [40], demonstrating a good description of different physical characteristics. 



IV. RELATIVISTIC CORRECTIONS AND THE WAVE FUNCTIONS AT THE ORIGIN 

With the use of the Eqs. (9) and (10) the kinetic energies cjj are calculated (see Table 11). Their numbers show 
that relativistic corrections are small in bottomonium: uji,{lS) — mi, ~ 180 MeV (~ 4%) and a bit larger, ~ 7%, 
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TABLE I: The vector constants As and Ajg^ (in MeV) (n/ = 3, 4, 5) for acrit = 0.605. 



rif :i 


i 




Ab 400 


372 


335 


Ams 271 


261 


245 



TABLE IL The kinetic energies ujq{lS){q = n,s,c) and ^(,(15') (in McV) for the static potential VB(r) (11) with n/ = 4, 
As = 372 MeV, acrit = 0.605. 

Meson B Bg Be 
mg 5 200 1410 

LOg - rUg 633 486 340 
rub 4790 4790 4820 
ojb - rub 42 46 110 



for the 25, 35 states. It is of interest to notice that a relativistic correction to the 6— quark mass in the Be meson 
is even smaller than in bottomonium, uib{lS) — rrn, — 110 MeV (~ 2%), while for a c— quark in Be the difference, 
a;c(l'S') — rUc — 340 MeV, is already ~ 25% (all numbers refer to UcHt = 0.605, n/ = 4). 
For the Bq mesons (g = n, s, c) relativistic corrections, WqilS) — rUg are given in Table II. 

In charmonium relativistic corrections to a c— quark mass are already ~ 13% for the ground state and ^ 17% for 
the 25 state. It is of interest to notice that they are even smaller, --^ 7% for a c quark in the D, Dg mesons (see Table 
III). 

For the analysis of HFS it is convenient to introduce the ratio ffs,, and also = g{bb),gc = g{cc), 

u^,{nS)o:,{nSy ^^^^ 



which directly enters the HFS (4): 



Ahf {nS) = - aHF (m) 9b, (nS) (19) 



and appears to be weakly dependent on small variations of the masses mi,m2, and other parameters of a static 
potential, which are compatible with a good description of the meson spectrum. Their values for the Bq mesons and 
bottomonium arc given in Tables IV and V. 

In general case the w.f. at the origin is sensitive to the choice of Asiuf), defining the OGE term, and may be ~ 1.5 
times larger, if the AF behavior is neglected [28]. In our calculations, when the AF property is taken into account 
and the same acrit{nf) is used for n/ = 3,4,5, the differences in the gB,,gD, values turn out to be small, < 5% for 
different nj, (see Table IV), with exception of gi, in bottomonium. 

Below wc give the values of gb{nf) for n/ = 3, 4, 5, which wore calculated with the same acrit — 0.605 . For this a^it 
and Mb — 0.95 GeV in two-loop approximation the constants Asirif = 3) = 0.40 GeV, ^sirif = 4) = 0.372 GeV, 
and Asiuf = 5) = 0.335 GeV are easily calculated (see Table I). The largest difference takes place in bottomonium, 
where 

gf,{nf = 3) = 0.213 GeV, gb{nf = 4) = 0.230 GeV, gb{nf = 5) = 0.258 GeV, (20) 



TABLE III: The kinetic energies ujq{nS){q = n,s) and lOc (in MeV) for the static potential Vb(»") (n/ =4) with the same 
parameters as in Table II. 



Meson 


D Ds 


cc(lS) 


cc{2S) 


mg 


5 200 


1410 


1410 


UJg — rUg 


542 400 


184 


245 


Wc — rric 


102 108 


184 


245 
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TABLE IV: The ratios t?s, (18)(in GeV) and \Ri{0)f (in GeV^) for the Bg{lS) mesons {acrit = 0.605, As(n/ = 4) = 372 MeV, 
As(n/ =3) =400 MeV). 

B Bs Be 
|i?(0)|2(n/ = 4) 0.510 0.558 1.742 
t,s^(n/=4) 0.165 0.168 0.212 
\R{0)f{nf = 3) 0.478 0.527 1.626 
gfl^(n/ = 3) 0. 158 0.161 0.198 



TABLE V: The ratios gb{nS) (in GeV), |i?„(0)|^ (in GeV^), and HPS AnpinS) (in MeV) for the 1S,2S,3S bottomonium 
states (n/ = 4,5, Ocrit = 0.605). 

IS 2S 3S 
\Rn{0)f{nf = 5) 6.476 3.398 2.682 
g(,(n/ = 5) 0.258 0.134 0.105 

AifF(n/=5) 71.1 36.9 28.9 
\Rn{0)f{nf = 4) 5.668 3.126 2.508 
g^{nf = 4) 0.230 0.127 0.100 

AHFjnf = 4) 63.4 35.0 27.6 



i.e. gb{nf) can change by ~ 20%. Therefore, if experimental HFS in bottomonium are known with a precision 

accuracy, then one can distinguish between cases with different nj. From the numbers given in (20) and taking 
UHF = 0.310, as for the Bg mesons, one obtains the fohowing HFS Ahf(&&) = M(T(9460)) - M{r]b{lS)): 

AHFirif = 3) = 58.7MeY; Anpinf = 4) ^ 63AMeY; AijF(n/ = 5) = 71.1 MeV (21) 

Notice that our value for the T(9460) — r]b{lS) mass difference with n/ = 3 has appeared to be in good agreement 
with the lattice calculations, also performed with n/ = 3, and where AHF{bb) = 61 ± 13 ± 4 MeV was obtained in 
Ref.[18]; 70 ± 11 MeV in [19], and a smaller splitting, AnFibb) = 54± 12 MeV, was calculated in [20]. (An accuracy 
of our calculations for HFS is estimated to be ±4 MeV (^ 5%), as it follows by varying different parameters, i.e. is 
better than in lattice calculations, where at present it is > 20%). 

For the Bg mesons their w.f. at the origin, the factors gs^, and the HFS arc given in Table IV for rij = 3,4. As 
seen from Table IV, for the Bg mesons a difference between gs^ {n / = 3) and gs^ {n j = 4) appear to be only by ~ 4%; 
also for a given nj a difference between gg and gg^ is small, < 2%. 

In bottomonium a difference between gb{nf = 3) and gb{nf = 4) is larger, ~ 12%, and in both cases they are 
smaller than gb{nf = 5); corresponding HFS are given in Table V. For excited states the T{nS) — r]b{nS) splittings 
{n = 2, 3) are given also for the coupling uhf = 0.310 and this choice seems to be a realistic, because in bottomonium 
a characteristic momentum weakly changes for excited states: 

LObilS) -mb = 185(5)McV, uj{2S) - mb = 195(5)McV, ^6(35*) - nib = 225(5)MeV. (22) 

In (22) theoretical errors given come from slightly different pole mass of a 6 quark for different n/. 

Finally in Tables VI, VII we give Qd^QDs for ground states, and also gdnS) = gcc{nS) for the IS" and 2S charmonium 
states. It is worthwhile to notice that the factor gcc{2S) is two times smaller than gccO-S), being one of the reasons 
why the mass difference -0(3686) — rjc{2S) is small. 

We would like to remind here that while the relativistic string Hamiltonian (7) is used, for a meson excited states 
may be considered on the same footings as a ground state till, if a single-channel approximation can be applied. 
Otherwise one has to use the multichannel Hamiltonian, also derived within FCM [41]. 



TABLE VI: The w.f. |i?(0)p (in GeV^) and the ratios go, go,, gc{lS), gc{2S) (in GeV) for the D, D, mesons and charmonium 
(in the static potential Vsir) cr = 0.18 GeV^, AB(n/ = 3) = 0.40 GeV, AB(n/ = 4) = 0.372 GeV, Mg = 0.95 GeV). 

D{1S) DsjlS) cc(15) cc{2SY 
\R{0)\^ 0.314 0.346 0.856 0.476 
g(nS) 0.379 0.379 0.340 0.174 



Tor the 2S charmonium state in the w.f. at the origin the S — D mixing between -0(3686) and i/>(3770) states is tciken into account 
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TABLE VII: The HPS (in MeV) in the Bq mesons with aHF{nf = 4) = 0.310, As(n/ = 4) = 0.372 GeV and aHF{nf = 3) = 
0.324, AB(n/ = 3) = 0.40 GeV. 

HPS B Bs Bc(lg) 

AHF{nf =4) 45.6 46.3 58.4 

AHF{nf = 3) 45.4 46.1 57.2 
AgF(exp) 45.78 ± 0.35 46.5 ± 1.25 abs 



TABLE VIII: The masses M{Bq{2S}}, M(B*(2S')), and the HPS M{B;{2S)-M{Bg{2S)), {q = n, s, c) (in MeV) for aHF{nf = 
4) = 0.310 and Ooritiuf = 4) = 0.605. 

Meson B(2S) B,(2S) B,{2S) 

gB^(2S) 0.124 0.1255 0.1353 

Ahf{Bq{2S)) 34.3 34.7 37.3 
M{Bq{2^So)) 5967 6040 6868 
M(B*(2^gi)) 6001 6075 6905 



V. HFS IN BOTTOMONIUM AND THE B^ MESONS 

Firstly, for bottomonium we compare calculated here HFS, Eq. (21), with experimental HFS, which follows from 
experimental mass of ?76(1S'): 

Ahf(66) = M(T(9460)) - M{r]b) = 69.9 ± 3.2 MeV [10]; 

Ahf (66) = 68.5 ±6.6 MeV [11]. (23) 

In the BaBar experiment the rjb mass is defined with a small errors, ±3.2 MeV, and even with a smaller error, 
< 1 MeV, the mass differences: M{B*) - M{B), M{B*) - M{Bs) are known now from experiment [7]. Such good 
experimental data allow to extract a coupling anp) with a good accuracy. For fitting procedure the important point 
is that in bottomonium a fitted value of anp cannot be larger that that for Bq mesons, since due to the AF property 
a coupling of a smaller size system (bottomonium) is typically smaller. For the Bg mesons and bottomonium the best 
fit is obtained for aHP = 0.310, later called "a universal" coupling. 

As seen from (21), in bottomonium full agreement with experiment is reached if Ohf = 0.310 and Uf — 5 arc 
taken in the static potential. It gives A.HF{bh) = 71.1 MeV, which is by 8 MeV larger than /\Hp{bb) = 63.4 MeV 
for n/ = 4. (We do not consider as unphysical a fit with n/ = 4 and aRF = 0.348 when agreement with experiment 
is also possible, since the value of this coupling is by 12% larger than that for the Bq mesons). Notice that the best 
description of the bottomonium spectra takes place also for n/ = 5 [40] . 

In Table V the HFS for the bottomonium 25', 35 states are given; they are equal 36(1) MeV and 28(1) MeV, 
respectively, being weakly dependent on n/ taken. 

The HFS of the Bq ground states are presented in Table VH, while in Table VHI the mass spittings for higher 
states, B{2S), Bs{2S), Bc{2S), are also given. 



For the B*{1S) — B{1S), B*{1S) — Bs{lS) splittings a good agreement with experiment is reached in two cases: 
with n/ = 4 and uhf = 0.310 (as in bottomonium) and with n/ = 3 and a bit larger uhf = 0.324 (see Table VII). 
First choice seems to be preferable as a universal one, but in any case a difference between two couplings is small, 
< 5%. Therefore one can speak about a imivcrsal uhf within 5% accuracy. To fix a preferable number nf for the 
B,Bg mesons one needs to use an additional information, like the decay constants etc. 

For the Be mesons our splitting, M{B*) - M(Bc) = 58.4 MeV for n/ = 4 and 57.2 MeV for n/ = 3, appears to be 
in agreement with the unquenched lattice calculations from [21], where the number 53 ± 7 McV is predicted. 

For excited Bc{2S) states our calculations give the centroid mass, Mcog{Bc{2S)) = 6893 MeV, and 
the B*{2S) - Bc{2S) splitting 37.3 MeV (see Table VIII), from which M{B*{2^Si)) = 6.902 MeV and 
M{Bc{2^Sq)) = 6865 McV. An accuracy of our calculations, performed in single-channel approximation, is 
estimated to be ±5 MeV , although for higher states an influence of open channel(s) may be important. 

From Tabic VIII one can see that for excited B,Bs,Bc mesons their HFS have close values, AHF{Bq{2S) ~ 
34—37 MeV. and the masses of singlet and triplet 25 states are also given in Table VIII. We would like to notice that 
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our HFS for the Bs racsons differ from those calculated in lattice QCD [22], where small Ahf{Bs{IS)) ^ 29.8±3.2 MeV 
was calculated for the ground IS* states (in our calculations it is equal 46.3 MeV), while on the contrary, in [22] for 
the excited 25* states a central value of the HFS, Ahf{Bs{^S)) = 56 ± 27 MeV, is larger than in our calculations, 
where this splitting is only 35 MeV. 

Our calculations have been performed in single-channel approximation with a string tension a = const = 0.18 
GeV^. For higher levels a influence of open channels can be taken into account, e.g. via a flattening of the static 
potential, and then the masses M{B{2S)),M{Bs{2S)) appear to be only by ~ 10 MeV smaller. The effect from open 
channels may be more important for the D{2S), Ds{2S) states and then a multichannel relativistic Hamiltonian from 
[41] may be used instead of the equation (7). 

Thus our analysis of HFS for the B{lS),Bs{lS) mesons and bottomonium shows that a good agreement with 
experiment is reached with a universal aHF = 0.310. This coupling is smaller than that in [2], [3] and corresponds 
to rather large renormalization scale, /i ~ 1.70 GeV. This value of the scale confirms existing interpretation of the 
spin-spin potential as dominantly a short-range perturbative one, thus justifying the use of the ^(r)-function. 

VI. LARGE HFS IN CHARMONIUM AND THE D, D, MESONS 

Experimental HFS for the D, Ds ground states are large, ~ 140 MeV, being three times larger than those for 
the B, Bs mesons. Let us firstly estimate these HFS, taking the factors go = go, = 0.379 GeV from Table VI 
[nj = 4), and the value anp = 0.310, as for the Bq mesons. Then for the ground states we obtain A.hf{D{1S)) = 
Ahf{Ds{1S)) = 104.4 MeV, which are by 35% smaller than experimental HFS [7]. 

In charmonium a discrepancy between calculated HFS with anF = 0.310, Ajjjr (15, cc) = 93.7 MeV,and experi- 
mental number is smaller, ^ 20%. If in charmonium first order correction (3), equal < 6%, is taken into account, 
then this discrepancy remains not small, ~ 15%. 

For the 2S charmonium states, a situation is different and full agreement with experimental HFS takes place, if a 
universal coupling anF = 0.31 is used. With the factor gc{2S) = 0.174 GeV from Table VI, one obtains 

Ahf{2S, cc) = M(3686) - M{r)c{2S)) = 47.9, MeV (24) 

coinciding with the experimental -0(3686) — rjc{2S) mass difference: A//i?(25, cc)|exp = 48 ± 4 MeV [7]. (To get this 
result we have taken into account the S — D mixing between tjj{3686) and ■ijj{377Q) with the mixing angle 9 ~ 11° 
[42]). Thus for the 25 charmonium states a universal coupling provides agreement with experiment. Therefore we 
expect that for the 25 states of other heavy-light mesons the coupling anp = 0.31 can be also used. 

However, it remains unclear what kind of corrections (~ 15% in charmonium ground states) have been lost in 
our analysis? (We remind that first order perturbative contribution gives only < 6%.) We assume here that such 
a contribution comes from nonperturbative spin-spin potential, and just due to nonperturbative correlators lattice 
calculations [18]-[20] give the mass difference, M{J/ijj) — M{ric{lS)), in good agreement with experimental number, 
equal 117 MeV. Such HFS, coming from nonperturbative spin-spin potential, can be taken into account also within 
FCM, where recently new results have been obtained for the vacuum correlation functions and correlation lengths 
[17], [43]. Detailed analysis of these effects will be considered in our next paper. 

However, we can use results from Ref.[25], thatfor higher states nonpcirturbtive contributions are much smaller than 
for the ground states; just for that reason we have obtained a good description of HFS for the charmonium 25 states 
with aHF = 0.31. Therefore we use here only perturbative part of HFS for higher states. 

It is of interest to notice that for D{2S) and Ds{2S) their w.f. at the origin and the quark kinetic energies can differ 
by ~ 10% , nevertheless the factors 5d(25) and gD,{2S) coincide within 4% accuracy, being equal 0.273 ±0.005 GeV, 
if linear confining term is taken in the static potential; it gives ~ 75 MeV for their HFS . If one takes into account 
flattening of confining potential, which is often important for higher states [44] , then a small decreasing of this factor 
takes place: gD{2S) ~ gD,{2S) = 0.264 ±0.004 GeV, where the error occurs due to possible different choice of n/ = 3 
or n/ = 4; in this case the HFS are equal 72 ± 3 MeV for the D{2S), Ds{2S) states. Calculated HFS are presented 
in Table IX. 

As seen from Table IX, the HFS of the D{2S), Ds{2S) mesons, ^ 70 MeV, are two times smaller than those for 
the ground states, i.e. for these states a picture is similar to that in charmonium, when the HFS for the 25 states is 
2.3 times smaller than the J/ip — rjcilS) mass difference. 

For the D{2S), Ds{2S) mesons our HFS appear to be significantly smaller than the Ds(2^5i) — £'s(2^5o) 
mass difference, predicted in [45], where it is equal 151 MeV, and even larger Ahf{D{2S)) = 188 MeV and 
Ahf{Ds{2S)) = 192 MeV were obtained in [46]. On the contrary, in our approach the HFS for the 25 states have 
appeared to be rather close to those from the GI paper [1], although in a static potential our and their sets of 
parameters are very much different, with an exception of the value of the string tension, equal 0.18 GeV^ in both 
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TABLE IX: The HFS (in MeV) of the D{2S), Ds{2S) mesons, and charmonium with aHF = 0.31; experimental HPS from [7]. 

Meson Ahf AHF{exp) 
D{2S) 72 ± 3 abs 
Ds{2S) 72 ±3 abs 
cc{lS) 93.7 116.6 ±1.2 
cc{2S) 47.9 49 ±4 



TABLE X: The masses of the singlet and triplet 2S states (in GeV) for D{2S) and Ds{2S) {aHF = 0.31). 



State this paper GI [1] RRS [45] MMS [46] 


D{2^So) 2,570 


2.58 abs 


2.483 


D{2^Si) 2.642 


2.64 abs 


2.671 


Dsi2^So) 2.664 


2.67 2.486 


2.563 


Ds{2^Si) 2.736 


2.73 2.637 


2.755 



calculations. Comparison our results and predictions from Refs. [1], [45], [46], where different relativistic models 
used, are presented in Table X. 

The masses, presented in Table X, were calculated with the use of different relativistic models, in particular the 
Spinless Salpeter Equation was exploited in [1], [45], and also in our calculations here, with the static potential 
described in Section 3. From Tabic X one can see that the singlet and triplet masses for the D{2S) and Ds{2S) 
mesons appear to be very close to each other in our calculations and in [1]. On the contrary, in [45] and [46] predicted 
HFS for the 25 states are 2.0 and 2.6 times larger than in our calculations, and due to this result, their masses of a 
singlet state, M{Ds{2^ So)), is by ^ 100 MeV and ^ 170 MeV, respectively, lower than in [1] and in our analyasis. 

For the triplet 25 states differences in predicted masses are not large and in [1], [46], and our calculations M{D* {2S)) 
and M{D*{2S)) lie in the range 2.64-2.67 GeV and 2.73-2.75 GeV, respectively. These predictions are in good 
agreement with the experimental mass, Mexp{Dl{2S)) = 2710 ± 2±7^ MeV [47], [48], while in [45] predicted mass is 
by ~ 70 MeV smaller. 

^From our analysis it follows that observation of the singlet states, D{2 ^5o) and Ds{2q), is crucially important for 
understanding of spin-spin interaction in systems of large sizes, in particular, it could clarify what is a characteristic 
value of the strong coupling in HF interaction. 

VII. CONCLUSIONS 

In our study we have used a conception of a universal HF interaction and observed that 

1. In the Bg mesons and bottomonium a good agreement with experimental HFS are reached if a universal coupling 
aHF = 0.310 is used in the HF potential. 

2. Just with the same coupling, uhf = 0.310, the -0(3686) — ric(2S) mass splitting appears to be in agreement with 
experiment. 

3. Calculated here mass splitting, M{B*) — M{Bc) = 57.9(6) MeV, gives the mass of unobserved yet B* meson, 
M{B*) = 6.334 ± 5 MeV. Our HFS is close to that in full QCD calculations, Ahf{Bc) = 53(7) MeV [21]. 

4. In bottomonium a full agreement with experimental mass of ?76(15) is reached only if in the static potential 
71/ = 5 is used, giving A//i?(66) = 71.1 MeV. For n/ = 4 and n/ = 3 calculated HFS are smaller, being equal 
63.4 MeV and 58.7 MeV, in agreement with the lattice results for n/ = 3. 

5. With aHF = 0.310 for the bottomonium 25 and 35 states the HFS, equal 36(1) MeV and 27(1) MeV, are 
predicted. 

6. The following masses of excited 5,(25) states are predicted: M{B{2S)) = 5967 MeV, M{B*{2S)) = 6001 MeV, 
M{Bs{2S)) = 6040 MeV, M{B;{2S)) = 6075 MeV, M{Be{2S)) = 6832 MeV, M(S*(25)) = 6870 MeV. 
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7. We predict that the mass differences for D*{2S) - D{2S), D*{2S) - 0^(23) are ~ 72(3) MeV, being smaller 
than in several other analysis. 

8. We expect that nonperturbative spin-spin potential gives not small contribution, ~ 15 — 30% to the mass 
splittings for J/V- - VcCi-S), D*{1S) - D{1S), and D*{1S) - D,{1S). 

For singlet states, D{2^Sn) and Ds{2^So), one cannot exclude that they have not small hadronic shifts due to 
coupling to open channels; if such hadronic shifts are small, then observation of their masses, close to the values 
^ 2.57(1) GeV and 2.66(1) GeV, can be considered as a crucial test of a universal character of spin-spin interaction. 
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VIII. APPENDIX A. THE POLE MASS OF A HEAVY QUARK 



The RSH (7) contains the pole mass of a heavy quark m2 = mq and of a hghter quark mi. For a light quark, m„ 
or rud, its mass is taken equal 5 MeV , while for a s quark mg = 200 MeV is used as in [26]. The pole masses of heavy 
quarks (c, 6) are defined as in pQCD, when the pole mass is expressed via the QCD current mass fhQ{fhQ), entering 
the QCD Lagrangian, and higher order corrections of the strong coupling [7]: 



niQipole) = mQ{mQ) 



l + ^^^+r. 
Sir 



V TT 



where the factor r2, 



r2(n/) = 13.4434- 1.0414 V (l - ll^) , 



(A.1) 



(A.2) 



depends on a number of flavors Uf through the sum, which goes from /c = 1 up to = n/ — 1. Due to this term the 
pole mass appears to be sHghtly different for different nj. Here in our calculations we take mf,(n/ = 5) = 4.823(3) GeV 
and mh{nf = 4) = 4.79 GeV; a small difference between them lies within a theoretical error, present in the current 
mass, TO(,(m(,) = 4.20 ± 0.07 GeV [7]. The mbipole) used here correspond to the value of the current mass frn,{fhb) = 
4.210 ± 0.015 GeV, which is within the conventional mmiber [7]. 

For a c quark the pole mass mc(n/ =4) = 1.41 GeV is used for charmonium, and the Be, Dg mesons; it corresponds 
to the conventional value, fhc = 1.24 ± 0.09 GeV [7]. For the D mesons mc{pole) = 1.39 GeV is used. 
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IX. APPENDIX B. THE QUARK SELF-ENERGY CONTRIBUTION TO A MESON MASS 

The nonperturbative quark self-energy (SE) contribution Ase to a meson mass M{nL) was calculated in FCM [32 
where the meson Green's function was defined in a gauge-invariant way. This correction, 



is negative, proportional to the string tension a, and a quark kinetic energy as w~ , which is rather large for a light 
quark and small for a heavy quark. In (A. 3) the factor r^q (a number) depends on a quark mass: = 1.0 for a light 
quark; rjg ~ 0.80 for a s quark; r]c — 0.40 for a c quark, and ~ 0.2 for a h quark [48], [32]. Notice that the SE term 
(A. 3) contains correct coefficient 1.5, instead of the coefficient 2.0 in [32]; the reasons for a change of this number is 
discussed in [49]. 

For a 6-quark the SE contribution is small (< (—1) MeV) and can be neglected; for a c quark its value is 

^ (—20) McV and it is convenient to include this small correction via a redefinition of a pole mass. For the ground 
states of heavy-light mesons the SE contribution, which comes from a light, is rather large, being ~ (—140) MeV, 
and a bit smaller, ~ —90) MeV, for a s quark. For higher states of heavy-light mesons the SE contributions are 
smaller, because for them the kinetic energy of a light (s quark), present in (A. 3), is larger. 



1.5a-?7(g) 



(A.3) 



